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1 Introduction 



Quantum chromo dynamics (QCD) is generally believed to be the correct theory for de- 
scribing the strong force between quarks and gluons. The property of asymptotic freedom 
makes it possible to use perturbation theory at large momentum transfer and allows for 
a precise description of many scattering experiments as carried out in the big accelerator 
facilities. 

As the fundamental theory for the interaction of quarks and gluons it also has to be 
able to describe bound states of quarks and anti-quarks, i.e. hadronic matter such as 
the proton, the neutron, pions etc. However, in this energy regime ordinary perturbation 
theory (expansion in the coupling constant) is bound to fail and QCD becomes extremely 
hard to solve. So far there are no promising analytical approaches at hand which would 
allow for a calculation of hadronic masses from QCD. It is the main goal of lattice gauge 
theory to numerically calculate hadron masses from first principles. 

Within lattice gauge theory QCD is not formulated in the continuum but on a discrete 
and finite space or space-time lattice. Hadronic masses can then be extracted from the 
decay of fermionic correlation functions. In Euclidean lattice gauge theory these correlation 
functions are given by Euclidean path integrals which in turn can be evaluated numerically 
by Monte Carlo methods. The fermionic degrees of freedom are formally integrated out and 
the Monte Carlo integration is carried out for a bosonic path integral only, which, however, 
contains the spectral determinant of the Dirac operator in the integration measure. 

The evaluation of fermionic determinants, which now has to be carried out for each 
update of the gauge field configuration, is computationally intensive, in particular for 
realistic, i.e. small, quark masses. Thus, in the past many studies have been performed 
in the so-called quenched approximation in which the fermionic determinant is neglected 
completely. This is equivalent to giving the quarks infinite mass or setting the number of 
flavours to zero. Large scale unquenched lattice calculations have only become available 
in recent years, and calculations with realistic quark masses will only be possible with the 
next generation of specialised super computers. 

Therefore, any other way of obtaining independent information on the spectrum of the 
Dirac operator, and thus the fermionic determinant, is of great interest. 

In the early nineties it turned out that chiral random matrix theory (RMT) describes 
spectral correlations of the QCD Dirac operator extremely well 0112] and can even predict 
the microscopic spectral density, i.e. the distribution of small eigenvalues of the Dirac 
operator, see j3] for an overview. However, the somewhat surprising information that the 
spectra of lattice QCD are, up to a certain scale, indistinguishable from the spectra of 
random matrices cannot be exploited directly in order to facilitate lattice calculations. 
The scale mentioned above is the equivalent of the Thouless energy in disordered systems. 
It was theoretically derived in jHEI and identified in lattice QCD data in [HIIZI- 

This situation is reminiscent of the situation in low dimensional quantum chaos. There 
short range spectral correlations of individual quantum systems can be described by RMT 
if the corresponding classical system is chaotic. In this context a two-fold role is played by 
semiclassical methods, in particular by the Gutzwiller trace formula |H]. On the one hand 
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they provide an explanation for the correspondence of classical chaoticity and quantum 
spectral correlations being described by RMT. On the other hand they also predict and 
describe deviations from RMT in long-range correlations, linking them to non-universal 
features of short periodic orbits jH]. Up to now, such a scale could not be identified in 
spectra of the QCD Dirac operator for a frozen, i.e. fixed, configuration of the gauge fields. 
We notice that the equivalent of the Thouless energy mentioned above is an effect due 
to the fluctuation of the gauge fields and can thus only be seen after averaging over all 
configurations, see the discussion in pUl lllj . 

Inspired by this analogy one may ask: Are semiclassical contributions the missing in- 
gredient which would make it possible to constructively use the RMT information when 
calculating fermionic determinants? As a first step towards an answer we develop semiclas- 
sical approaches to the Euclidean QCD Dirac operator and in particular discuss qualitative 
features of the classical dynamics arising in this context. Notice that the word "semiclas- 
sical" in this context always refers to asymptotic statements about the spectrum of the 
Dirac operator - technically a problem in single particle quantum mechanics rather than 
in quantum field theory - which is not the same as loop expansions which are also called 
"semiclassical" in quantum field theory. Thus, our approach is in a similar spirit as works 
relating the spectral analysis of the QCD Dirac operator to the theory of disordered systems 

laEj. 

This article is organised as follows. In section El we review some basic formulae and 
discuss the semiclassical structure of the QCD Dirac operator. In section IHl we briefly 
sketch a strategy for deriving trace formulae which we will follow in the subsequent sections. 
The discussion of semiclassical approaches to the Dirac operator in Abelian gauge fields 
presented in section|3]serves as prerequisite for our semiclassical analysis for the QCD Dirac 
operator which follows in section |3| The latter contains the main results of this article 
identifying three semiclassical parameters and discussing the classical dynamics arising in 
different (combined) semiclassical limits. Some details left open in sections |3] and El are 
solved by our study of the squared Dirac operator in section IHl In section [71 we discuss 
whether and how our theory can be used for describing features of Dirac spectra close to 
zero virtuality. Section |H1 illustrates our theory for an explicit example. We conclude with a 
discussion of our findings and by indicating possible future directions of research in section 

El 



2 The QCD Dirac operator 

The free Euclidean Dirac operator describing massless spin 1/2-particles reads 

D = ^l,d,. (2.1) 

We adopt the summation convention over repeated Greek indices from 1 to d, the number 
of space-time dimensions. The 7-matrices satisfy 

{7^, 7.} = 25^,. (2.2) 
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Describing massive particles simply amounts to adding — im to ()2.ip . Also note that in 
the context of lattice gauge theory often the anti-Hermitean operator iD is called Dirac 
operator. In dimension c? = 4 we will later explicitly use the chiral representation, 



where 1„ denotes the n x n unit matrix and cr is the three-vector of Pauli matrices, 

a: = (? J) > (° 7) , ^3 = (J . (2.4) 

In this representation 75 = 71727374 reads 

When now introducing a non-Abelian gauge field we put special emphasis on the ap- 
pearance of h and its consequences for semiclassics. In this way we will identify the most 
natural asymptotic treatment from the perspective of semiclassical physics. However, by 
introducing fields with an /i-dependent magnitude alternative options are also possible 
and we will remark on those in the appropriate places. Moreover, we will explain which 
situations in standard QCD language correspond to the scenarios discussed. 

A non-Abelian gauge field Afj_[x) is introduced by minimal coupling, 

^ = 7/. (y^/. - hgA^i^)^ . (2.6) 

Notice the appearance of h, together with the coupling constant g which turns the covariant 
derivative into 

D^ = d^- igA^ix) , (2.7) 

D = —ik'-f^Dfj^. This is different from the Abelian case, i.e. quantum electro dynamics 
(QED), where the minimal coupling prescription reads 

^d, ^ ^d^-eA,ix). (2.8) 

Here we denote the coupling constant, i.e. the electric charge of the fermion, by e (we set 
c = 1) and thus the covariant derivative reads 

D^ = d^-i^A^. (2.9) 

The reason for the different appearance of h in these two cases is that non-Abelian fields 
couple to themselves. More precisely, when writing down the QCD-Lagrangean which 
upon variation yields both, the Dirac equation and the classical Yang-Mills equations for 
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A^, the latter would contain a self-interaction term which would explicitly depend on h if 
the covariant derivative ()2.7p had the same /i-dependence as ()2.9p . Since this cannot be 
true for a classical equation, formula ()2.7|1 is the correct choice for non-Abelian fields. To 
illustrate this consider the field strength tensor deriving from ()2.7j) . 

= d^A, - d,A^ - ig[A^, A,] (2.10) 

which does not contain h. On the other hand an Abelian field does not couple to itself, the 
last term in ()2.10|1 vanishes, and therefore in QED a covariant derivative of the form ()2.9|1 
is allowed, because it does not lead to an ^-dependence of the field strength. The very same 
mechanism is responsible for the well-known fact that one can have elementary particles 
with different electric charges but that all particles which couple to the non-Abelian colour 
field do so with the same coupling constant, i.e. they all have the same colour charge. 

Notice that the observation described between ()2.6jl and here holds true as long as all 
/i-dependencies in the formulae are displayed explicitly, i.e. quantities such as g, A^ or F^^ 
do not depend on h. In particular, the argumentation as laid out above is independent 
of the scaling properties of the QCD-action. If one, e.g., rescales the fields according to 

= A'^/V^ then eqs. (EH) and (ITTn|) will read D = i^, {jd^, - h^'^gA'^) and F^^ = 
y/hdf^A'^ — ^/hdiyA'^ — ihg[A'^, Al], respectively. Thus, we have formally produced powers 
of h in unfamiliar places. However, as long as the original fields A^ are of order 1 then 
the rescaled fields A'^ are of order h~^^'^, and therefore the ^-dependence of the couplings 
between fermion and gauge field and of the gauge field to itself are as before. 

The situation changes if we, instead of just rescaling the fields, consider fields whose 
order of magnitude is /i-dependent. For instance, a gauge field of order 1/h gives rise to a 
Dirac operator in which the coupling of fermion and colour field has the same /i-dependence 
as in QED. In QCD such a strong field is called an external colour field. 

In the situation, however, which was described between ()2.6p and ()2.10j) . the electro- 
magnetic fields are external fields whereas the colour fields are microscopic or dynamical 
fields. Since the /i-dependence of the latter is chosen such that the classical field equations 
are /i-independent it is also common in QCD to speak of a "classical gauge field" in this 
context. From the point of view of high energy physics it may appear slightly inconsistent 
to discuss external electro-magnetic fields and microscopic colour fields in the same con- 
text. From the point of view of semi classical physics, however, we have treated both types 
of fields on the same footing. 

In the following we will concentrate on the situation with microscopic gauge fields, i.e. 
on the Dirac operator (j2.6p where all /i-dependencies are displayed explicitly. A treatment 
of external colour fields would lead to different semiclassical asymptotics. 

In order to shed some more light on the physics behind the /i-dependence discussed 
above, decompose the non-Abelian field in terms of the generators X°- of the gauge group 
G, say SU(A^), 

A,= ^-A';,X\ (2.11) 

where summation over the repeated Latin index a is from 1 to the dimension of the Lie al- 
gebra. The X" are traceless, Hermitean x A^-matrices satisfying the Lie algebra relations 
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[X^^X''] = r^'^X^, (2.12) 
with structure constants /°'"^, and are normalised according to 

tT{X'^X'') = 26ab. (2.13) 

If we now view 

:= ^X" (2.14) 

as the quantum observable describing the colour degrees of freedom (of the fermion) the 
Dirac operator ()2.6p takes the form 

D = ^,(p,-gC''A;{x)) , (2.15) 

with the colour and momentum operators and p^, representing the quantisation of some 
classical observables and p^. 

The point of view adopted in (j2.14j) is typical for internal, i.e. microscopic, degrees 
of freedom, a familiar example being the non-relativistic spin operator s = ha/ 2, which 
has the same structure as ()2.14j) with X'^ replaced by the Pauli matrices, the generators of 
SU(2). 

For later reference let us also introduce the matrix valued function on classical phase 
space, 

D{p. x) = 7, (p, - ^X^A'^ix)^ , (2.16) 

from which the QCD Dirac operator can be obtained by replacing with —itid^. In 
microlocal analysis or Wigner-Weyl calculus ()2.16|) is referred to as the Weyl symbol of the 
Dirac operator (j2.6|) . which in turn can be obtained from its symbol by Weyl quantisation, 

{m){x) = -i-j / [ d(p, ^) e^r>.i^^-y^) vi;(^) d^^d^^. (2.17) 

Wigner-Weyl calculus is a particularly useful tool when studying semiclassical asymptotics. 
In a setting where the semiclassical limit is identified with h ^ one would classify the 
terms in ()2.16|) according to their /i-dependence as the principal symbol 

Do{p,x)=^^p^ (2.18) 

and the sub-principal symbol 

D^ip,x) = -^^,X^A^^ix), (2.19) 

respectively. Eventually we will also use the notation 

sjmh[D]{p,x) = D{p,x) (2.20) 
for the Weyl symbol of an operator. 
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3 Semiclassical trace formulae 



Before we go into details about the semiclassics for the QCD Dirac operator let us say a 
few words about semiclassical trace formulae in general and briefly sketch one method for 
deriving them; for details, however, we refer to the cited literature. 

We are interested in the spectrum of the Hermitean operator D. For simplicity assume 
that the spectrum is pure point, i.e. we have a set of eigenvalues and a complete 
ortho-normal set of corresponding eigenstates 

b^r^ = \n^n. (3.1) 

Our main focus lies on the spectral density 

p(A) = 5^5(A-A.), (3.2) 

n 

which is gauge invariant. In order to derive a semiclassical expression for p(A) consider the 
evolution equation 

d 

ih—^{x,t) = Dm{x,t). (3.3) 

Note that the time parameter t is not the physical time but an auxiliary variable. The 
physical time is already included in the components of x and we are dealing with the 
Euclidean Dirac operator. Likewise the spectral parameter A is not an energy but referred 
to as virtuality. Now define the evolution kernel K{x, y, t) by 



^{x,t)= / i^(a;,2/,t)vl/(y,0)d^, (3.4) 
which has the spectral representation 

K{x, y,t) = Y, ^n(x) n(2/) e"^'"* . (3.5) 

n 

Obviously, K{x,y,t) also has to solve ()3.3j) with initial condition 

K{x,y,0)=6{x-y). (3.6) 

By Fourier transforming the evolution kernel and taking the trace on L^(M'^) (g) C"' ® C^, 

where J denotes the dimension of the representation of the gauge group, we obtain the 
spectral density, 

tr ^ / [ K{x,x,t)e^^'dtd'^x = p{X) (3.7) 

Here tr denotes the trace over the matrix degrees of freedom. 

In oder to obtain a semiclassical approximation for the spectral density, one can begin 
with a WKB-type ansatz for the time evolution kernel. 



K{x,y,t) = j^ I [ao(x,e,t) + /iai(x,e,t) + ...]et(^(^'«'*)-«^)d'^e- (3.8) 
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Inserting into ()3.3p and sorting by powers of h one finds a sequence of equations which 
can be solved order by order yielding S,aQ,ai, .... In leading order one always finds a 
Hamilton- Jacobi equation for the phase S, 

A{V,S,x) + — = 0, (3.9) 

with a classical Hamiltonian A given by an eigenvalue of the principal symbol of D. Clas- 
sical Hamilton-Jacobi theory now tells us that the solution S of ()3.9j) generates classical 
dynamics from the phase space point V^S) to (VxS, x) in time t, showing that the in- 
tegration parameter ^ of the ansatz ()3.8j) plays the role of an initial momentum for the 
classical system. 

In order to derive a semiclassical approximation to p(A), one also needs to determine 
the leading order amplitude oq which is fixed by the next-to-leading order equation. The 
latter, usually referred to as transport equation, has the following structure, 

\dt ^ ""y ° 2\dXf,dp^ dp^dpu dXf^dxu J ^ "' 

The reader easily verifies this structure by explictly doing the calculation for an operator 
of his choice. A derivation of the general result can, e.g., be found in appendix E of [T^ . 

The first bracket on the l.h.s. of (j3.10p is a derivative along the trajectory gener- 
ated by the A-dynamics, whereas the second term, roughly speaking, measures the be- 
haviour of neighbouring phase space points. Without additional terms ()3.10jl is solved by 
^^detW^Jdxd^, see, e.g., appendix B of [12] for a compact derivation. If, besides the terms 
displayed explicitly in ()3.1()|1 . further contributions show up in the transport equation then 
they represent the transport of internal degrees of freedom (such as spin or colour as we 
will see below) along the trajectory of the flow with Hamiltonian A. 

Therefore, in order to find the full classical system corresponding to the quantum 
Hamiltonian D one has to (i) determine the Hamiltonian(s) A and (ii) carefully analyse all 
contributions to the transport equation. 

Having determined S{x,^,t) and a(x, ^,t), i.e. having obtained a semiclassical ap- 
proximation to the kernel K[x,y,t), a trace formula can be derived in a straight-forward 
manner by inserting everything into ()3.7p and evaluating all integrals in leading order with 
the method of stationary phase. 

The stationarity conditions for the x- and ^-integrals read 

V^5(x,^,t) = ^ and V^S{x,^,t) = x . (3.11) 

According to classical Hamilton-Jacobi theory this means that both initial and final mo- 
mentum as well as initial and final position of the trajectory generated by S have to 
be identical. Thus, only phase space points that lie on periodic orbits contribute to the 
semiclassical expression for p(A). 

A special role is played by the periodic points of period zero which are given by the 
whole hypersurface of constant virtuality A, 

nx:={ip,x)\A{p,x) = \}. (3.12) 
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Since their action is also zero they yield the only non-oscillating contribution to the spectral 
density and thus constitute the mean density, often called Weyl term, 

with a possible multiplicity factor deriving from the internal degrees of freedom such as 
spin or colour. 

Finally we obtain the following general structure for a semiclassical trace formula, 

p(A) ~p(A) + ^^^(A)e^^-W. (3.14) 

7 

Here 7 labels both, isolated periodic orbits and larger families of periodic points, like, e.g., 
Liouville- Arnold tori in integrable systems. The amplitudes ^^(A) are derived by keeping 
track of all contributions in the various stationary phase approximations involved. 

If one is interested in the precise mathematical meaning of this distributional identity 
and in an absolutely convergent version of the trace formula, which can be used for nu- 
merical calculations, it is convenient to multiply the expressions with a test function in t 
before taking the Fourier transforms, see e.g. [T!^ ITU ITK] 



4 Semiclassical parameters in the Abelian case 



In this section we discuss semiclassical approximations to the Dirac operator in Abelian 
gauge fields. We will keep the presentation short and closely follow similar studies for 
the Dirac Hamiltonian which were carried out in [THl [T3] . however, pointing out small 
differences which are due to the fact that we are dealing with the Euclidean Dirac operator 
instead. The results obtained here will also be needed for our discussion of the non- Abelian 
case in the following section. 

Consider the equation of motion ()3.3p for the time evolution kernel. 



where the derivatives in the Dirac Hamiltonian, 



^ = 7m ( - eA^{x] 



(4.1) 



(4.2) 



with Abelian A^, act on the first argument of K. Inserting an ansatz of type ()3.8j) with 
scalar phase S and matrix-valued amplitudes into the evolution equation we find in 
leading orders 



— + D{V,S,x) 
ot 



— + D{y^s,x) 



d 



Oq = 



«i + ( ^ + 7m'9m ) Oo = 



(4.3) 
(4.4) 
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where 

D{p,x)=( °. ^4-i^7r\ 

^ \^7r4 + 10-77 y ^ ^ 

is the (principal) symbol of D and vr^ '■= Pfi — denotes the kinetic momenta. For 

()4.3|) to have non-trivial solutions the term in square brackets must have an eigenvalue 
zero. The eigenvalues of D{p,x) are given by 

A±(p,x) = ±y^=:±A, (4.6) 

both having multiplicity two. We collect the corresponding eigenvectors columnwise in the 
2 X 4-matrices 



I f lo \ .... 1 r 



7r4 — ICTTT 



n(p.x) = ^^^j. V4p,.) = ^^[_\J. (4,7) 
Thus, the solvability condition for ()4.3jl yields the Hamilton- Jacobi equations 

A±(V.^±,x) + ^ = 0, (4.8) 

and the general solution of ()4.ip is a superposition of terms with positive (+) and negative 
(— ) virtuality. 

Equation (j4.8p alone does not solve (j4.3|) . but in addition the leading order amplitude 
has to satisfy Da^ = A^a^. This is guaranteed by the following ansatz, 

4{x, t) = V^iV.S^, x) b^{x, e, t) Vli^, y) (4.9) 

where the 2 x 2- matrices h± have to be determined by (j4.4p . To this end we multiply (|4.4|) 
with Vl.iVxS'^^x) from the left and V±{^,y) from the right, yielding 

Vl{V^S^,x) (^7A + ^) ^±(V.^±,x) 6± = , (4.10) 



smce 

Vl{V,S^,x) — + D{V,S,x) 
I at 

After a lengthy calculation, which is sketched in appendix one finds 



(4.11) 



Vil4^ + Vh,V^(d,b^) = (g^ + S„j 6± =■ b± , (4.12) 
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where the dot in ()4.12|) denotes a derivative along the trajectory generated by S^. In 
addition we have introduced the electric and magnetic components, E and B, of the field 
strength F^j,, according to 



F 








B3 


-B2 


-Ei\ 









Bi 


-E2 




B2 


-Bi 





-E-s 


V 


El 


E2 


E3 


/ 



(4.14) 



Since we already know how to solve a transport equation of type (I3.10|) the product ansatz 



det T^-TT d± 



dxd^ 



with a 2 X 2 matrix d lends itself to simplify the transport equation to 



d. 



± 



le 
2A± 



(T{E±B)d± = 0. 



(4.15) 



(4.16) 



This equation describes the transport of the spin degrees of freedom along the trajectory 
determined by the Hamiltonian A^. Obviously d± takes values in SU(2) and in the trace 
formula the contribution of each periodic orbit is weighted with the trace of the corre- 
sponding rf-t, i.e. with a character. 

Equation ()4.1f)|l can be mapped from SU(2) to 5^ by looking at the time evolution 



of the expectation value s of the spin operator in an arbitrary state u G 
s± = u'^ dl.^(Td±u - as induced by ()4.16|) . 



I.e. 



s± = s±x —{E±B) 



(4.17) 



This equation describes classical spin precession, i.e. it is a Euclidean analogue of the 
Thomas- or Bargman-Michel-Telegdi(BMT)-equation ^71^]. Although ()4.17p looks like 
an equation for the three-vector s classical spin precession actually takes place on the 
sphere S"^ since total spin, i.e. is conserved. The two-sphere in turn is a symplectic 
manifold and ()4.17|1 defines a volume-preserving fiow on it. These facts together justify 
the notion of "classical spin dynamics" in this context. Mathematically speaking, we map 
the equation from (the representation of) the group to its coadjoint orbit, see e.g. [T^OU] . 

Had we dealt with a particle with higher spin from the beginning we would have ob- 
tained a similar spin transport equation as (j4.16p with only cr replaced by generators of 
a higher dimensional representation of SU(2) and d± now taking values in that represen- 
tation. The weight factor tTd± in the trace formula would still be a character and the 
analogous mapping to the sphere would lead to exactly the same classical spin precession 
equation. The character entering the trace formula is always completely determined by 
classical spin precession [T^ . 
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We have thus identified the total classical dynamics arising from a semiclassical analysis 
of the Euclidean Dirac operator as a combination of the Hamiltonian fiows with Hamilto- 
nians A^[p,x) accompanied by spin precession along the orbits. Since there is no back- 
reaction of spin dynamics on the Hamiltonian part the total dynamics can be formulated 
as a skew product fiow, either on R^d x SU(2) or R^'^ x S"^ flUU^. 

So far we have discussed what happens in the single semiclassical limit h ^ 0. However, 
also the limit of large spin can be considered as a semiclassical limit, cf. the so-called kicked 
top [22] • If one simultaneously lets h —>■ and s — oo, where 2s -|- 1 denotes the dimension 
of the representation of SU(2), such that the product hs is kept non-zero and finite, also 
the back-reaction of spin on the translational degrees of freedom shows up in the classical 
picture, see e.g. [211123 1211 • We emphasise that, although claimed otherwise in [23], even 
for Hamiltonians linear in the spin degrees of freedom semiclassical asymptotics can only 
then display both, spin dynamics and back reaction, simultaneously in leading semiclassical 
order if one considers the large spin limit. This says, however, nothing about the possible 
practical use of this type of approximation even when the actual value of s is rather small. 
The situation here is somehow reminiscent of the limit of large colour, N oo, [2ZI, which 
yields valuable insights although we are mostly interested in = 3. 

We can almost write down the Hamiltonians A^{p,x,s) for the combined dynamics 
already with the information gathered so far. Omitting the spin-dependent terms it has to 
reduce to A'^{p,x), i.e., formally, A='^(p, a;,0) = A=''(p, x), and it has to give rise to the spin 
precession fl4.17|l . The relativistic Pauli Hamiltonian, 

A^,,i,(p,x, s) = A±(p,x) - ^^^^s{E{x) ± B{x)) , (4.18) 

fulfils these requirements, but so does, e.g., the alternative square root type Hamiltonian 

A±rt(p,x,s) = ±y/ {p^ - eA^{x)){p^ - eA^{x)) - 2es{E{x) ± B{x)) (4.19) 

Both types of Hamiltonians agree in the limit of small spin contribution (for illustration 
one may formally consider the limit s ^ 0) but in general they lead to different back reac- 
tions of spin on the translational degrees of freedom. This difference becomes particularly 
important for small virtualities A. In section [HI we will dicuss a simple method for deciding 
which Hamiltonians to use, without explicitly developing a full symbol calculus for the 
combined limits. 

5 Semiclassical parameters in the non-Abelian case 

With a discussion of semiclassical parameters and limits of the QCD Dirac operator and 
the different classical dynamics arising in this context the present section contains the 
central results of this work. We perform our analysis along the same lines as laid out in 
sections [2] and [3] and build on the results obtained in section [3J 

Consider the Dirac operator (j2.6j) with Weyl symbol (j2.16|) . The time evolution kernel 
K{x, y, t) is now a 4 J x 4 J matrix, where J denotes the dimension of the representation 
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of the gauge group, i.e. J = 3 for QCD with SU(3) gauge fields in the fundamental 
representation. Inserting an ansatz of type ()3.8|) into the equation of motion ()4.H) with 
Dirac operator ()2.6j) yields in leading orders 



OS 
'dt 



ao = - 



dS 
'dt 



ai 



d_ 
di 



7A + ^i(Vx5,a;) ao = 



(5.1) 
(5.2) 



where we have used the notation for the principal and sub-principal symbol which was 
introduced in eqs. ()2.18|) and ()2.19|) . The principal symbol has eigenvalues 



A±(p, x) = ±y%pjl =: ±A 
with corresponding eigenvectors collected in the matrices 



po+ia-p j ) 
A / 



1 / pa-icrp 



V2\ 



A 
-l2 



(5.3) 



(5.4) 



As in section 01 eq. ()5.H) demands that the phase of the ansatz ()3.8j) solves a Hamilton- 
Jacobi equation, 

A±(V.^±,x) + ^ = 0, (5.5) 
at 

and suggests the following ansatz for the leading order amplitude, 

4{x, e, t) = y±(V,^±, x) 6±(x, t) Vii^, y) . 
The projected transport equation for the 2J x 2J matrix h± reads 



(5.6) 



(5.7) 



which, using eqs. fj4.12|) and ()4.13|) with the substitution vr^ p^, cf. (|5.4p . simplifies to 
1 / 92A± 92A± d^S^ 



2 \dx^dp^ dp^dpudx^dxu 



(5.8) 



As in the Abelian case we separate the translational part according to (j4.15|l and obtain 
the following equation for the 2J x 2J matrix d^-, 



2 dpf, ^ 



(5.9) 



In contrast to the Abelian case this equation does not involve the spin but the colour 
degrees of freedom. Accordingly we will refer to it as colour transport equation. 
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As in the case of spin transport we obtain classical dynamics by looking at the equation 
of motion satisfied by the expectation value C° of (7" = 

c"^ = -|^r'^4^% (5.10) 

which we call colour precession. Equation (j5.1(Jj) is the colour part of the Wong equations 
[2H] to which we will come back later. 

As in the case of spin precession discussed in the preceeding section there are certain 
conditions restricting the possible values which the variables C"" can assume, thus confining 
the dynamics ()5.10p to a compact manifold: The (representations of the) Casimir operators 
of the gauge algebra are constants of motion for and from those derive constants of 
motion for the precession equation ()5.1()|1 . In the case of SU(2) there is only the quadratic 
Casimir operator (total spin in the preceeding section) which confines the precession to the 
sphere S"^ = SU(2)/f/(l). If the gauge group is SU(3) then we have two Casimir operators, 
one quadratic and one cubic in the generators (or components of classical colour). The 
dynamics of the is thus reduced from (8 generators) to a six dimensional manifold, 
the flag manifold = u(fy|7^(Y)' ^-S- 122] • For the gauge group SU(A^) we would 
have = SU(A^)/U(1)'^~^, instead. In all cases these are maximal coadjoint orbits [TU] . 
which are not only even dimensional but naturally endowed with a symplectic structure, 
thus constituting the classical phase space for internal degrees of freedom such as spin or 
colour. 

Having understood ()5.9p as transport equation for the colour degrees of freedom and 
characterised the underlying classical phase space and dynamics we can now ask ourselves 
why the spin degrees of freedom do not show up at this level of the semiclassical treatment, 
neither in the Hamiltonians nor in the transport equation. The answer is that spin 
and translational degrees of freedom are coupled by the gauge fields and thus only via 
the internal colour degrees of freedom. With both, spin and colour, being internal degrees 
of freedom a coupling between them, which has to involve the product of s = | cr and 
C"" = ^X"-, is automatically at least of order h^. Therefore, it does not enter the leading 
order phases and amplitudes of ^ — > asymptotics. 

Comparing with the results of the preceeding section we should expect a spin precession 
equation like ()4.17|) with the electric and magnetic fields replaced by their non-Abelian 
analogues. At this point we can thus guess the following semiclassical hierarchy (which we 
will confirm to be correct in the following section): In pure ^ — > asymptotics the phase 
of semiclassical approximations is determined by free translational dynamics alone. The 
leading order amplitude is affected by the colour transport along particle orbits, whereas 
there is no back-reaction of colour onto the translational degrees of freedom. Spin shows 
up only as an ^-correction to the amplitude. While spin precession is driven by both 
translational and colour dynamics it does not act back on either of them. Thus, we have 
a double skew product structure. 

Back reaction can be forced to show up explicitly in the semiclassical approximations 
by considering combined limits. To this end choose a J dimensional unitary irreducible 
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representation of the gauge group and consider the combined hmits h and J ^ oo 
with the product hJ kept constant. This will lead to colour entering on the same level as 
the translational degrees of freedom, i.e. we will have to deal with the minimally coupled 
classical Hamiltonians 

A±(p,x,C) = ±^{p, - gA-{x)C-){p, - gA'>^{x)C^) . (5.11) 

The coupled classical dynamics arising from these Hamiltonians are known as the Wong 
equations ^Hl ISO] . Spin will enter on the level of the transport equation as for pure ^ — 
asymptotics in the Abelian case. Thus we have moved from a double skew product structure 
to an ordinary skew product. We mention in passing that by performing the additional 
limit J ^ oo with hJ fixed we have changed the order of magnitude of the term ghA^, 
appearing in the Dirac operator ()2.fij) . from h to 1. Thus, this situation may be physically 
related to that of an external colour field in the language of QCD, cf. the discussion 
following ()2.10|) . Mathematically, however, the scenario introduced here is different. 

If we go even one step further by taking the triple limit ^ — 0, J — oo and s ^ oo 
with the products hJ and hs fixed we will find fully coupled Hamiltonian dynamics on the 
total phase space M^'^ x x S"^ (for SU(3)-gauge fields). As for the relevant Hamiltonians 
we have to solve the same problem as at the end of section |3] With the knowledge obtained 
so far it could be either of Pauli type ()4.18|) or of square root type ()4.19|) . with A^, E and 
B replaced by their non-Abelian analogues. 

6 The squared Dirac operator: 

Confirming the semiclassical hierarchy 

The following study of the squared Dirac operator serves two purposes. On the one hand 
we prove that the semiclassical hierarchy conjectured in the preceeding section is correct 
and on the other hand we determine the functional form of the classical Hamiltonians 
corresponding to the Dirac operator in simultaneous semiclassical limits. 

We calculate the square of the Dirac operator ()2.6p and determine the Weyl symbol of 
for different symbol calculi. The principle underlying this approach is that in a symbol 
calculus there exists always a so-called Moyal product which expresses the symbol of the 
product of two operators as an asymptotic expansion in the semiclassical parameter(s) in 
terms of the symbols of the individual operators. The leading term in this expansion, i.e. 
the principle symbol of the product, is always given by the product of the principal symbols. 
Thus, from the appearance (or absence) of certain dynamical variables at given order in 
the symbol of the squared Dirac operator we can conclude at which order these variables 
appear in the symbol of the operator itself. Moreover, the principal symbol of the squared 
operator allows us to draw conclusions about the functional form of (the eigenvalues of) 
the principal symbol of the operator itself. 

The square of D is most conveniently calculated by decomposing the products 
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and D^Dy into their symmetric and antisymmetric parts, 

= -h'^^^^D^D, = ({7^, 7,} + [7^, 7,]) {{D^, D,} + [D^, DJ\) 

= {{l,.lu}{D,,D,} + Yl,,^,][D,,D,]) (6-1) 



/lU j 1 



where on the last line we have used ()2.2|) and the definition ()2.10|) of the non-Abelian field 
strength tensor, F^j^ = ^[D^, D^]. With the representation (j2.3p we have the commutators 

[74, 7] = 2i and [7,, 7;^] = 2i£,,, °j , j, fc, / = 1, 2, 3 . (6.2) 

Introducing colour-electric and colour-magnetic fields as in ()4.14|1 we finally obtain 

= -reD.D^ - teg f '^^^ + ^) \ . (6.3) 



(t{B-E)J ■ 

Thus, the matrix-valued Weyl-symbol of reads 

symh[b^]{p, x) = (p, - hgA,){p, - hgA,) - h'g (^'^^^ + . (6.4) 

From this we can easily read off at which order in h the different degrees of freedom will 
enter a semiclassical approximation. To this end recall that A^ = ^A'^X"', E = ^E"'X"' and 

B = ^B'^X"-. At order only the translational degrees of freedom show up in symb[-D^]. 
The colour degrees of freedom, X"-, appear for the first time at order h^, whereas the spin 
degrees of freedom are absent unless we proceed up to order h^. The first two observations 
are in agreement with the semiclassical analysis of the non-Abelian Dirac operator, and 
the last one provides the missing element in order to prove the semiclassical hierarchy 
anticipated at the end of sectional 

If, instead of Wigner-Weyl calculus for the translational degrees of freedom only, we 
used a symbol calculus which also maps the internal matrix degrees of freedom, spin and 
colour, to classical variables, i.e. symb[|cr] = s and symb[|X"] = C"", then the symbol of 

D'^ reads 

syrah[D']ip,x) = (p, - gA^C^P, - gA^C) - 2g + - E^)C') " 

(6.5) 

In the simultaneous limit h 0, J ^ 00 and s — > 00 with hJ and hs fixed this total 
symbol consists of a (diagonal 2x2) principal symbol only, i.e. there are no higher order 
terms in any of the three semiclassical parameters. The linearity of (j6.5|) in the spin degrees 
of freedom shows that the eigenvalues of the symbol of the Dirac operator D itself have to 
be of square root type ()4.18p rather than of Pauli type ()4.19p . 
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We remark that in the Abehan case ()6.4|) reads 



symb[D^](p, x) = {p^ - eA^){p^ - eA^) - he 



(t{B + E) 







a{B - E) 



) 



(6.6) 



which is consistent with spin appearing in the leading order transport equation in pure 
h ^ Q asymptotics. 



We have motivated this study with the question whether and how semiclassics can be of 
use for the understanding of spectral properties of the QCD Dirac operator. Of interest 
are here in particular spectral functions averaged over an ensemble of gauge fields, as they 
appear in the calculation of path integrals in (lattice) quantum field theory. An important 
example is the averaged spectral density, 



where the action S[A\, and thus the integration measure, can, e.g., be just the Yang-Mills 
action (quenched approximation) or the bosonic part of the full QCD action, including 
fermionic determinants. 

Prominent features of (p(A)) are the so-called chiral condensate, a non-zero value at 
virtuality A = 0, and a universal functional form for small values of the virtuality. Due to 
the Banks-Casher relation jSI], see also [321 El El], (p(0)) is proportional to the expectation 
value {ipip) of the quantised quark fields '0 in the ground state. Thus, a non-zero value 
indicates the spontaneous breaking of chiral symmetry. Moreover, after suitably rescaling 
the virtuality with the chiral condensate, the microscopic density becomes universal and 
can be calculated in chiral RMT. Notice that the chiral condensate is not determined by 
the number of exact zero modes, which is a topological invariant of the Dirac operator, 
but it arises, in a suitable limit, from the accumulation of small but non-zero eigenvalues. 

In this section we investigate whether and to what extent the different semiclassical 
approaches characterised in section El are able to explain the formation of a chiral con- 
densate on the level of the Weyl term (jH.lHj) . which is the semiclassical description of the 
mean spectral density. For these considerations we restrict the x-integration to a subset 
V C M'^ with finite volume V and since the density is symmetric about zero it is sufficient 
to consider only positive A. 



In pure ^ — ^ asymptotics the periodic orbit structure of the trace formula and the hy- 
persurface Vtx fl3.12|) determining the Weyl term are derived from the translational degrees 



7 Mean density in stochastic fields 




(7.1) 



7.1 h 
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of freedom only. Moreover the translational dynamics are extremely simple, namely free. 
Spin and colour enter only as multiplicity pre-factors, 



J{2s 



r- [ [ s{,/pijr,-x)d''pd''x 



^ 27r'^/V(2g + 1)V 
(2vr^)^r(f) 

Thus, the Weyl term reduces to the free mean spectral density and is in particular inde- 
pendent of the gauge field configuration. For h = 1, dimension d = 4, J = N = 3, the 
number of colours, and spin s = | the mean density reads 

m = ^,x'- (7.3) 

7.2 /i ^ 0, J ^ oo 

If we take this combined limit then colour precession ()5.10|) appears on the same level as 
translational dynamics. The Hamiltonians are not only functions of p and x but also of 
the classical colour degrees of freedom C. In order to be able to integrate over this larger 
phase space we need a parameterisation of the colour part. If the gauge group is SU(A^) (in 
a faithful representation) then C has A^^ — 1 components. Colour dynamics, however, lives 
on the A^(A^ — l)-dimensional manifold F^. Let C{^) be a parametrisation of F^, then the 
correctly normalised integration measure, which complements d'^pd'^x/{2'n'hY stemming 
from the translational degrees of freedom, is 

d^(^-i)e (7.4) 



|F^i 

with |F^| = Jpjv (i^(^-i)^ being the volume of the flag manifold. Thus the Weyl term now 
reads 

m = JwJ^^P^[ [ [ 5iA^p,x^CiO)-\)d'pd'xd^^^~'k. (7.5) 



|F^| {2TThy 



However, with A"*" = {Pn — 9^1i{x)C"'){p^ — gA^^{x)C^), see (jSHIl, after a simple shift of 

variables in the p-integrals for fixed x and C, Pfi ^ Pfi + 5'Aj^(x)C", integration over the 
colour degrees of freedom becomes trivial and ()7.5p reduces to ()7.2j) and thus once more to 
the free result. 
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7.3 h ^ 0, J ^ oo, s ^ oo 

In this triple limit all degrees of freedom - translational, colour and spin - appear on the 
same footing in the Hamiltonian, 



= ^{P, - 9A-{x)C-){p, - gA>l^{x)C') - 2gs{E\x) + B\x))C- . (7.6) 

In order to calculate the Weyl term we also need to parametrise the phase space of spin, 
5^, which we do in spherical coordinates denoting the solid angle by uj. The correctly 
normalised measure is ^^d^cj, with volume 15"^! = 47r. The mean density thus reads 



(2s+ 1)J 



47r|F^|(27r;i) 
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P(A)= . "Zm\:S.., / / / / 5(A+(p,a:,C(0,sM)-A)d'^pd'^xd^(^-i)ed^^ 



(7.7) 

As before we can shift the integration variable in the p-integrals in order to remove the 
explicit appearance of the gauge potentials AJ^. However, through the field strengths E 
and B the expression still depends on the gauge fields and the integration over the internal 
degrees of freedom does not become trivial. Yet we are able to calculate the mean density 
if we consider the average over an ensemble of gauge fields, which is the function we are 
interested in anyway, 

(p(A))= ! -p{\)VA. (7.8) 



Using the following property, 

5(A+-A) = 2A5(A+'-A2), (7.9) 
and employing the Fourier representation of the (5-function we have to calculate 

(p(A)) = -%tll:^ / [Iff / e^(^^'-^')M'^pd'^a;d^(^-i)ed^^dtDA. 

(7.10) 

For simplicity we calculate this expression using stochastic fields. More precisely, we assume 
locally independent Gaussian fluctuations with the same variance a for all components of 
E and B. According to a relation derived in appendix El for Weyl terms this is equivalent 
to averaging over constant random fields, i.e. 

...VA^ ^4mxr^ / • • • e-(^''^"+^"^")/(2'^') ^^NiN-i)j^^zNiN-i)j^ 

(7.11) 

Now the total exponent is quadratic in E and B and an average over the fields yields. 



1 



^-i2gs(E''+B"')C"'t g-(£<'i;°+B°B<')/(2cr2) ^ZN{N-1) ^m{N-l) 



(27ra2)3A^(A^-i) ^7^^2) 
= exp {-Aa'^g^s^C'^CH'') . 



19 



Since and C^C^ correspond to the quadratic Casimir operators of SU(2) and SU(A^), 
respectively, they are constants, i.e. they depend only on s and J but not on uo and ^. 
Hence, 



An'^/^VJ{2s + 1) A 



TTV^ Jo 



where we have introduced the abbreviation v := \^ ag\s\y/C"-C°-. This parameter, being 
proportional to the variance a and the coupling constant g, is a measure for the strength 
of the fields. The field free situation corresponds to = and one easily confirms that 
in this case (I7.13j) reduces to (j7.2p . Similarly, for large A the integral expression grows 
proportional to A"^"^ and we once more obtain the free density ()7.2|) . For arbitrary A the 
remaining integral in ()7.13|) can be expressed in terms of generalised Laguerre functions. 
The most important observation is that for A ^ the integral converges to a constant. 



e pd-i = S — Z_ r(f). (7.14) 



1 r -p' H-^ . (2t;2)'^/4-i/2 
Jo 

Thus, for small virtualities the mean density now grows linearly instead of being propor- 
tional to A*^"^ as in the previous cases. For instance, for d = 4, where the behaviour changes 
from cubic to linear, this means a dramatic increase in the number of small eigenvalues. 
If the triple limit discussed here was related to the scenario of a strong external field then 
the linear density for small A could be interpreted as the density within the first Landau 
band, cf. [3^. In any case, the behaviour of the density resulting from our semiclassical 
approach shows some remarkable features which we discuss in the following. 



7.4 Discussion of Weyl terms 

We have calculated the (averaged) Weyl terms for the spectral density of the QCD Dirac 
operator in all 3 different semiclassical limits introduced in section |3| In two cases the 
resulting leading order mean density is just the mean density for the free Dirac operator. 
Only in the triple limit, — > 0, J ^ oo and s ^ oo, have we observed a dependence on 
gauge fields and internal degrees of freedom. In particular we have derived an increase in 
the number of small eigenvalues. 

In this last case we have replaced the QCD or Yang-Mills action in the path integral 
over the colour fields by a Gaussian measure, thus neglecting details of the gauge dynamics. 
The success of random matrix models and in particular related work on stochastic field 
theories fOl makes us believe that our results are nevertheless relevant for QCD. 
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To further study the accumulation of small eigenvalues as borne out by ()7.13p we have 
to examine the integral expression 



$,(A):=-^/ e-^(^'-^')'/-Mp. (7.15) 
y2'KV^ Jo 



Its value in dimension c? = 4 is given by 

$4(A) = / e-^(^'-^')' dp 




(7.16) 



For large A the integral is dominated by the last term, i.e. it grows like A^/2, which restores 
the A^-behaviour of the free mean density ()7.3|) . On the other hand for small A the integral 
is determined by the first term, giving rise to a linear spectral density, i.e. (for d = A, 
J = 3, s = I and h = 1) 

It is interesting to see how this expression changes under the variation of external param- 
eters. If one wants to consider finite temperatures in a field theoretical setting one has to 
choose an asymmetric subset V C M^, say a box with lengths with fixed L4 -C Li, L2, L3. 
Then the inverse of L4 is essentially the temperature, see e.g. [33]. Having one smaller 
dimension sets a natural scale for the small eigenvalues, namely \ k, 2tt / (cf. the eigen- 
values of the free Dirac operator in a box with lengths L4 <^ Li, L2, L^). Since it is this 
accumulation of small but non-zero eigenvalues which we want to investigate further, it is 
instructive to look at the averaged spectral density on this scale. 



Lj/ (27r)3/2 L4 V 2t^2^^ 
In fig. [T]we plot the scaled density 



31/V27r 



(P) (7.19) 



function of the scaled variance 



observing a curve which is reminiscent of a critical phenomenon with the spectral density 
itself playing the role of the order parameter. The behaviour for C > Cc would be interpreted 
as an indication for a non-zero density for small A, which, on the other hand, vanishes 
(exponentially) for ( < (c- Thus, we are tempted to view Q as a value indicating a phase 
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transition. In terms of the original quantities, the variance measuring the (couphng) 
strength of the gauge fields, and the inverse temperature L4, this implies the following. 
For fixed temperature, on the one hand, the phase transition would occur at a critical 
strength of the gauge fields with a vanishing order parameter for weak fields. For fixed f , 
on the other hand, we would observe the phase transition for a critical temperature with 
a non-vanishing order parameter at low temperatures only. 

We find it remarkable that our simple semi- 
classical argument is capable of showing a be- 
haviour which seems to hint at a critical phe- 
nomenon. In view of this we are tempted to 
put forward the following speculation. Accord- 
ing to the Bank-Casher relation [21] a non-zero 
averaged density at A = results from the for- 
mation of a chiral condensate. However, in the 
derivation of the Banks-Casher relation the limit 
A ^ may only be considered after one has first 
taken the infinite volume limit V 00, for a 
suitably normalised expression, and then the chi- 
ral limit m — >■ of vanishing sea quark mass(es). 
If one interchanges the latter two limits, i.e. if one 
performs m —* before V ^ 00 then the chiral 
condensate vanishes, see e.g. [3^1 • In our semi- 
classical calculation we do not have a mass parameter, which we could vary accordingly. 

Nevertheless, we find an accumulation of small eigenvalues hinting at a critical phe- 
nomenon. Could it be that this relates to the formation of a chiral condensate? - However, 
we do not want to conceal that our discussion only takes into account the leading order 
Weyl term. Higher order semiclassical corrections to the mean density may contribute 
where the leading order term vanishes, as we will demonstrate for an example in section 
ISJ and the above discussion also ignores the periodic orbit contributions to the spectral 
density. Recall that a trace formula provides a decomposition of the density of states p into 
a mean term p and a periodic orbit sum posc, cf. flH.14j) . of which only the latter oscillates 
as a function of the virtuality A. It is now natural to expect that after averaging over the 
gauge fields only the non-oscillating Weyl-term contributes to (p). The above discussion 
is based on this tacit assumption. In general we have 




Figure 1: Graph of the scaled density 
r(C) = Cexp(- 1/2(2), gqg_ ([7191) and 
(|7.2U|) . which we speculate to be a measure 
for the chiral condensate. 



(p(A)) ~ (p(A)) + (po.c(A)) 



(7.21) 



and it is not guaranteed that (Qosc) = holds in a short virtuality interval close to A = 0. 

If besides the formation of the chiral condensate one also wants to explain the universal 
microscopic density, characteristic of the chiral ensembles of random matrix theory and 
observed in lattice calculations, cf. [3], a theory involving (correlations in) the periodic 
orbit contributions will be required. The analysis in j^ZlEHlEni! where the emergence of 
universal microscopic densities is discussed within a graph model should be viewed as a 
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guideline which, combined with our semiclassical approximations, would put the semiclassi- 
cal understanding of spectral correlations in QCD on a similar level as in (low dimensional) 
quantum chaos. 



8 Example: Fermions in SU(2)-fields on 

As an illustration for the structure of semiclassical trace formulae and for the calculation 
of some of the contributions we discuss the example of the Dirac operator on a two- 
dimensional torus with constant SU(2) gauge fields. With "constant field" we actually 
mean constant potentials yl^, which, due to the non-Abelian character, can give rise to 
a non- vanishing field strength, cf. ()2.1()j) . For this scenario we can analytically calculate 
the eigenvalues and derive an exact trace formula. The contributions to this trace formula 
are then compared to the corresponding semiclassical expressions. In order to keep the 
presentation simple we only discuss the case with fixed representations for the internal 
degrees of freedom, i.e. pure h asymptotics. 

Consider a two-dimensional Euclidean Dirac operator, in external SU(2) fields, 

D=(%-^9A,a]^,. (8.1) 



.i 2 

In two dimensions the 7-matrices can be chosen of type 2x2, e.g 



7.= I? I). 7.= (? -'|. (8.2) 



JO/' Vi 

More precisely, we should write 

7^ = I2 ® a'^ (8.3) 

and replace cr in ()8.1|) by cr ® I2. As configuration space we choose a two dimensional box 
with lengths and periodic boundary conditions, 

*(X1 +Li,X2) = ^(Xi,X2), ^{XuX2 + L2) = ^ixi,X2) , (8.4) 

i.e. we put the system on a torus T^. With the ansatz \E^(x) = u exp{j^p^x^) the boundary 
conditions require 

27r^ 

Pfj^ = , n^^l (no summation convention!) (8-5) 

and the Dirac operator reduces to an ordinary 4x4 matrix, which has to be diagonalised. 
The eigenvalues are most conveniently determined via the square of the Dirac operator. 
By a calculation similar to that in section IHl one finds 



= \Ip,P^^ + ^9'A,A^ ± y h^g^p,p,A^A, + ^g\A, x A^)^ (8.6) 
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with as in ()8.5p . The spectrum is once more symmetric about A = and we only show 
the positive eigenvalues. 

For the following we concentrate on the special case with 



(8.7) 

in which yet none of the four contributions in ()8.fij) vanishes. Introducing the abbreviation 

a := hgA (8.8) 

we have 



, / „ / a? , a 



A„ = y p^.p^. + y ± Y + X " V ^^^^ ^ T ^ 2 ■ 

The spectral density (for positive virtuality) thus reads 

p{X)=Y. [5(A-A+) + 5(A-AJ] . (8.10) 

When the spectrum is already known exactly a trace formula can usually be derived by 
employing the Poisson summation formula which expresses a sum over integers by a sum 
over the Fourier transformed addends, 

E/H = E/ f{n)e'^"'^^^d''n. (8.11) 

Doing this for the spectral density ()8.1()j) and changing variables from to = L^n^j {2nh) 
(no summation convention) we have 

P(A) = 7^$^/ KA-A+)+5(A-A-)]ei^^^"dV, (8.12) 

where we have introduced 

Qfi '■= kf^Lfj_ (no summation convention) . (8.13) 



Introducing radial coordinates, p = A/pf+~p|, q = Y/AjfZf~+"AJ|Z|, we obtain 



with 



\Hp)--=^p' + ^±^ (8.15) 
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The (5-functions select p = y^A^X^^oy and p = A(A + a) in the first and second term, 
respectively, and the 0-integral yields a Bessel function. Hence, 



P(A)-|^[e(A-a)(A-i) + (A + i 



[e(A - a) (a - I) Jo (I v/A(A^) + (a + ^) Jo (| V^AIaT 



(8.16) 



where we have separated the mean density, 

^w-i§h^-°'(^-i)+(^-i)]- ("^) 

which derives from ki = k2 = 0. Notice that - since a is of order h, cf. (|8.8|) - in leading 
semiclassical order the mean density is given by 

P(A)~^A (8.18) 

which agrees with (j7.2p with d = 2, J = 2 and, formally, s = 0, since there is no dynamical 
spin in 1 + 1 dimensions. For small A, however, the exact mean density reads 

^<^) = (^ + 2) • 

giving rise to a non-zero value at A = which cannot be seen by leading order semiclassical 
asymptotics. 

Before we can compare the periodic orbit sum with semiclassical theories we have to 
expand the result asymptotically for h —>■ 0. To this end recall that a is of order h. Using 
the asymptotic behaviour of the Bessel function we obtain 



Jo (I V'A(A±^) ~ J cos (|v/A(AT^- ^) 

V^"^ / y 7rg^A(A±a) '4/ 

2h (qgA\ /q 71 
cos COS —A 



(8.20) 



nqX \ 2 J \h 
and thus the semiclassical periodic orbit sum reads 

Let us analyse some contributions. The rapidly oscillating term, the last cosine, contains 
the argument qX/h. The geometric length of a periodic orbit on the torus with winding 
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numbers ki and k2 is given by g = ^ykfLf + k^L"^ and with the Hamiltonian ()5.3|) the 
action of this orbit reads 

jp^Ax^ = \q. (8.22) 

The colour field shows up only in the argument of the other cosine, which has to derive 
from colour precession. For SU(2) fields the colour transport equation ()5.9|) reads 

d-i^A^(Td = Q. (8.23) 
2A 

With the choice ()8.7j) this becomes 

aA 

d-i^p^a''d = (8.24) 

which we have to integrate with initial condition d{0) = I2 up to the period of a periodic 
orbit on the torus. Due to the Hamiltonian ()5.3|) the period equals the geometric length 
of the orbit and thus the solution is given by 

d = exp f-i^p^a^^q] . (8.25) 



The trace of this expression, 

, , ^ ( gA^jpi^^ \ (gAq\ 

trd = 2 cos I ^ q 1=2 cos I I , (8.26) 

enters as a weight factor in the periodic orbit sum, and indeed gives rise to the aforesaid 
cosine factor. The remaining factors can be calculated with standard methods, see e.g. 
|4m HT| for the general case or |12[ section 3.6.2] for a related example. 



9 Conclusions and outlook 

We have discussed the semiclassical structure of the QCD Dirac operator, and in partic- 
ular the interplay of three semiclassical parameters, namely Planck's constant h, and the 
spin and colour quantum numbers s and J, respectively. This situation allows for various 
semiclassical scenarios, with combined semiclassical asymptotics considered. We have en- 
countered a rich family of classical dynamics of translational, colour and spin degrees of 
freedom, whose mutual coupling depends on how the semiclassical limit is taken. 

The infiuence of these different types of dynamics in semiclassical trace formulae has 
been discussed and, in particular, we have analysed the behaviour of the Weyl term, the 
mean density of states, in different semiclassical scenarios. Based on this analysis we have 
critically evaluated which of the semiclassical scenarios has the potential of describing 
the spectrum of the QCD Dirac operator near zero virtuality, leading us to a speculative 
discussion of the mechanism behind the chiral phase transition. We certainly do not want 
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to overstate this speculation, a definite statement requires further work, as indicated in 
section I7.41 There are various directions of research which would naturally continue the 
present analysis. 

So far we have mainly discussed the Weyl term, which in a trace formula gives rise to 
the mean density of states, but not the periodic orbit sum which is responsible for spectral 
correlations. An analysis based on periodic orbits should, e.g., lead to a semiclassical theory 
for the universal microscopic spectral density of the QCD Dirac operator as described 
by chiral RMT. Moreover, such an approach would also describe deviations from RMT 
behaviour on large spectral scales, cf. saturation effects as described in jH], and thus 
potentially provide the missing link asked for in the introduction, which would make it 
possible to directly use RMT information when calculating fermionic determinants. Here, 
one should keep in mind that the equivalent of the Thouless energy j3J Ej sets another 
scale. It is not present in Dirac spectra for frozen gauge fields, but it is an ensemble effect 
resulting from the propagation of the gauge fields. 

Besides the various semiclassical scenarios which we have described in this article there 
is an additional strategy for taking the semiclassical limit of multi-component wave equa- 
tions. In this approach one does not treat the matrix degrees of freedom, i.e. colour and 
spin, dynamically but rather considers polarised Hamiltonians, which describe a particle 
with the spin or colour projection locked to the "direction" of the external gauge field, see 
e.g. |1Z1 UHl Uni mi- Such an approach may also prove useful in the case of the QCD Dirac 
operator. 

In lattice gauge theory, which we have referred to in various places, the Dirac operator 
is implemented as a difference instead of a differential operator. This has consequences 
which could also be analysed within the semiclassical picture. On the one hand discreti- 
sation leads to a modified dispersion relation, i.e. to different classical Hamiltonians A. 
Roughly speaking, the momenta are replaced by suitably normalised sines of momenta, 
which in the semiclassical picture changes both, Weyl terms and the periodic orbit struc- 
ture. On the other hand the discretised theory lives in a finite dimensional Hilbert space. 
In a semiclassical context the dimension of this Hilbert space also becomes a semiclassical 
parameter (cf. the theory of quantised maps, jl^l; see also [IH] for an overview), which 
would allow for the continuum limit to be discussed on a semiclassical footing. 
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A Projected transport equations 

In order to calculate the projected transport transport equations ()4.1()|1 for the 2x2 
matrices b± we have to evaluate the expressions, 



db± 



Vl^V± and Vh^d^Vl 



(A.l) 



We begin with the terms where the derivatives act on h±: 



7r4 — ICTTT X + ICTTT 

(flo + ica) + (ao — icra) 



A 



A 



1 



2A 



(2ao7r4 + 2a7r) 
9A 



A 



(A.2) 



Due to normalisation, Vl.V± = I2, the second contribution is trivial, and together with the 
first one yields 



db+ 
'dt 



Vh,V+{d,b, 



(A.3) 



where the dot denotes a derivative along the flow with Hamiltonian A. The remaining two 
terms yield 



dt 



V2 



A2 



[714 — ICTTT ) 







x^""^ dt 



2A2 



diz . 97r4 

lVr4Cr— ICTTT—— 

at at 



(A.4) 



and 



1 
2 
1 
2 



(So - io-V 



7r4 + iCTTT 



-— 7r^(ao7r^ - icrV7r^)(7r4 + Icttt) + ^(^q - icrV)(7r4 + i(nz) 



(A.5) 
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Using the Hamilton- Jacobi equation ()4.8|) we can derive the useful relation 

dS_ 



which we now use "backwards" , 

1 



A3 



1 / (97r4 d7z 



Hence, 



Is? \^~dt + TT— + ICr I — X TT 

1 

+ A 



dt dt ' " dt V dt 

(9o7r4 + icr^oTT — icrV7r4 + Vtt + icr(V x tt)) 



and added to ()A.4|1 a couple of terms drop out, 



+ 2A^ ~ + V X tt) 



2 [(?3:^t?p/j c)pfj,dpi^ dx^dxi^] 



With 



SoTT - V7r4 = 9o(V5 - gA) - V{doS - gA^) = -g{doA - VAq) = -gE 



and 



V X TT = V X (VS - gA) = -gV x A = -gB 
the non-scalar terms finally are given by 



2A 



-a(B + E) 



The analogous calculation for 6_ reads 



ao - icra \ ( — 



7r4 + ICTTT 



A 



[ao — icra) — [ao + icra) 



-I2 

7T4 — iCTTT 

A 



-^(2ao7r4 + 2a7r) 



A 
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where the dot now denotes a derivative along the flow with Hamiltonian A . The other 
two terms yield 



d_ 

dt 



Vl-V^ = ^ ^ ^-^'^ 



( ~dt '""^ _ 7r4-igr7r Stt^ 
A A3 "a* at 



— (7r4 + icr7rj 



9t 



icr 



dt 



2A2 



17r4Cr— h ICTTT— h ICr 77 X — ) 



dt 



dt 



dt 



(A.14) 



dtJ 



and 



V2 



(9o + io-V)^^ 



(9o + i(TV) 



7r4 — ICTTT 

A 



-^7r/,((9o7r^ + i(7V7r^)(7r4 - i(77r) + -^(Sq + icrV)(7r4 - iaiz) 



Again we use the Hamilton- Jacobi equation (j4.8p . 

as 

"Tt 



-dT = di^^^ = = -^''^ =5.A=X^.vr., 



(A.15) 



(A.16) 



concluding that 



vh.d^vl 



1 

"A2 



dni . dn^^ . diT diz 



9^ 



X TT 



+ — (9o7r4 — icr^oTT + icrV7r4 + Vtt + icr(V x tt)) 



(A.17) 



and together with ()A.14j) we obtain 



vl^v^ + vh,d,vl 



2h? dt 2A 



icr 

2A 



(^oTT — V7r4 — V X tt) 



(A.18) 



With (jA.lOjl and (jA.lljl the non-scalar terms in this case read 

-^(t(B -E) = ^^criE - B) . 

2A ^ ^ 2A- ^ ' 



(A.19) 
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B Local Gaussian fluctuations vs. constant random 
fields 

We show that averaging Weyl terms over stochastic fields - more precisely, independent 
locally Gaussian fields - is computationally equivalent to averaging over constant Gaussian 
fields. 

When averaging a mean density p, given by a Weyl term, over fields, say B{x), it is 
crucial that the Weyl term is given by an integral over the position variable x, see (j3.13|) . 
Therefore, let us now consider expressions of the form 



/:= j J f{B{x),x)dxVB {B.l] 



with a Gaussian measure VB. Think of the functional integral as defined by a suitably 
normalised continuum limit oo of its discretised analogue, 

with lattice points xj and B{xj) =: Bj. Then we obtain 

i.^Y. ^1^"^^ - / ^i^dB, (B,3) 

J=l V 
In the continuum limit we have thus derived the relation 




f{B{x),x)dxVB = j J f{B,x)—^=^dxdB, (B.4) 

i.e. when interested in local Gaussian fiuctuations in section 17. 3| we may average over 
constant random fields instead. 
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